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^ Abstract 

The General Rotorcraft Aero mechanical Stability Pro- 
gram (GRASP) was developed to analyse the steady-state 
and Unearned dynamic behavior of rotorcraft in hovering 
and axial flight conditions. Because of the nature of prob- 
lems GRASP was created to solve, the geometrically non- 
linear behavior of beams is one area in which the program 
must perform well in order to be of any value. Numerical 
results obtained from GRASP are compared to both static 
and dynamic experimental data obtained for a cantilever 
beam undergoing large displacements and rotations caused 
by deformation. The correlation is excellent in all cases. 

Nomenclature A' 

l ] 

b£ = basis vectors at the deformed beam tip j 

bf = basis vectors at the beam root « 

c = beam width ! 

Cii - direction cosine matrix relating the beam root 
and tip 

e = error in beam cross-sectional measurement 

E = modulus of elasticity 

7a = geometrical cross-sectional property 

7« = geometrical cross-sectional property 

p = projection of bf onto the bf - bf \t 
t — beam thickness 

U — strain energy 

a = material nonlinearity coefficient 
0 — torsional deflection which was experimentally 

measured 

k = curvature of the beam 

4>i = Rodrigues parameters 

$ = load angle 


Introduction 

The General Rotorcraft Aeromechanical Stability Pro- 
gram (GRASP) 1 is capable of treating the nonlinear static 
and linearized dynamic behavior of structures represented 
by collections of rigid-body and beam elements that may 
be connected in an arbitrary fashion and are permitted to 
have large relative motions. GRASP was developed primar- 
ily for analysis of the steady-state and linearized dynamic 
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behavior of rotorcraft in the hovering and axial flow flight 
conditions. 

GRASP combines multibody and finite-element tech- 
nology by taking the strong points from each area and in- 
tegrated them together into a single, comprehensive pack- 
age. GRASP differs from standard multibody programs by 
considering flexible-body and aeroelastic effects, including 
simple, nonlinear, unsteady aerodynamics. GRASP differs 
from standard finite-element programs by allowing multi- 
ple levels of substructures in which the substructures can 
move and/or rotate relative to others with no small-angle 
approximations. This capability facilitates the modeling 
of rotorcraft structures, including the rotating/nonrotating 
interface, and details of the blade/root kinematics for var- 
ious rotor types. An overview of the features of GRASP 
as a program can be found in Ref. I, while details of the 
analysis are given in Ref. 2. The theoretical basis of the 
analysis are addressed in companion papers. 3,4 

Because of the nature of the problems GRASP wis cre- 
ated to solve, the geometrically nonlinear behavior of beams 
is one area in which the program must perform well in or- 
der to be of value. The main structural element in GRASP 
is the aeroelastic beam, a geometrically nonlinear beam el- 
ement based on the kinematics, internal and inertial forces 
of Ref. 5 and the aerodynamics of Ref. 2. The element 
kinematics are valid for small strains and large rotations, 
but shear strains are neglected. The beam element degrees 
of freedom include (as generalized coordinates), rigid-body 
translations and rotations thst are associated with the stan- 
dard cubic and linear shape functions. These account for 
the usual twelve degrees of freedom found in beam finite- 
element analyses. In GRASP, however, an arbitrary num- 
ber of additional generalized coordinates associated with 
higher-order polynomials® are also included, allowing the 
element to have a number of degrees of freedom that is ar- 
bitrarily more than the usual twelve, following the so-called 
*p-versioa" of the finite element method. 7 * 1 

This paper is to present numerical results from GRASP 
for comparison with static and dynamic experimental data 
for large deflections of an end-loaded cantilevered beam. 
(The experimental apparatus and procedures are only 
briefly described in this paper. Details of the experiment 
may be found in Ref. 9 and 10.) Determination of the 
beam properties for input into GRASP is described in de- 
tail. GRASP result* are then presented along with results 
from previous analyses which are shown for comparison. 
While the present results do not exercise many of the fea- 
tures and power of GRASP, they do serve to validate much 
of the code dealing with the beam element 1 * ability f x> model 
highly nonlinear behavior. 
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Experiment 

An experiment done it Princeton University®* 10 (un- 
der Aaroflightdynamics Directorate sponsorship), was se- 
lected as a test case with which to validate GRASP. This 
experiment consisted of measuring the static deformation 
and fundamental flatwise and edgewise natural frequencies 
of a uniform, nonrotating, cantilever beam with a mass 
4 *toched to the tip (Fig. 1). The beam was slender and 
sufficiently flexible to undergo large displacements (still at 
small strains) due to the preeence of the tip mass. Beam 
load angle and mass of ths tip weight were varied through- 
out appropriate ranges. 

Ths beam wac instrumented with strain gages mounted 
at the root in the proper orientation to measure flatwise 
and edgewise natural frequencies. The end of the beam 
was securely mounted in a precision indexing chuck that 
provided a stable mount and accurate, repeatable angular 
settings. The static vertical and horixontal measurements 
were made with a caliper scale, measuring from a flat table 
with a reference grid affixed. The static torsional measure- 
ments were made with the aid of lightweight reference rods 
attached along the length of the beam and perpendicular 
to It. (See Ref. 10 for complete details.) 

The beam was fabricated from 7075 aluminum. The 
mass density was assumed to be 2.626 x 10- 4 Ib.-see. 3 /in.\ 
Ths length of the beam was measured to be 19.985 in. 
The t hickn e ss and width of the beam were measured at 
0.1251 in. and 0.4999 in., respectively. Assuming a gravita- 
tional constant equal to 380.089 in./sec. 3 , the mass per unit 
length was determined to be 1.6424 xlO"* Ib.-sec . 3 /in. 3 . 
The mase moments of Inertia were 2.1420 xlO”* Ib.-sec 3 
(flatwise) and 3.4204 xl0~ 7 Ib.-sec. 3 (edgewise). 

Determination of the appropriate values of bending 
stiffness proved to be more difficult. Both static and dy- 
namic predicted behavior are sensitive to the value of the 
stiffnesses, therefore stiffnesses were determined as accu- 
rately as possible. Attempted Inference of equivalent beam 
properties from classical linear formulas for deflection vs 
load for ths two uncoupled cases, load angles of 0 and 
9CT , yields contradictory information — even when only 
small deflections are considered. At a load angle of 0* (the 
edgewise- bending case), linear theory is too stiff. But, at 
a load angle of 90* (the flatwise-bending case), linear the- 
ory is too soft. Similar contradictory information results 
when attempting to correlate natural frequencies of the un- 
loaded beam according to linear-beam theory with exper- 
iment. This suggests that there is no one value of E that 
will yield accurate flatwise and edgewise bending stiffnesses 
if the measured cro ts - ea ction dimensions are taken as exact 
and the theory is assumed to be linear. With the failure 
of linear theory, we turned to a simple, nonlinear, planar 
elasttca model. 11 Even then, a standard value of the modu- 
lus of elasticity and measured croes-sectional dimensions in 
a standard elastka model, only fair agreement is obtained 
with planar, experimental, static deflections. 

With the aforementioned problems in mind, a more in- 
novative approach was called for. First, it was determined 
that only static data should be used because of supposed 
accuracy. Second, because transferee deflections for the un- 
coupled, no-tip- mass cases were recorded to be sero (even 


though non- sero deflections should have been measureable), 
it was necseeary to assume that the experimental data have’ 
deflections for the notip- mass case subtracted out. This 
latter point is not explicitly stated in Ref. 9 or 10, but 
from the present investigation, appears to be true. Finally, 
to *Uow for other effects not present in the simple elastic* 
model, an extra parameter in the form of a material nonlin- 
earity coefficient a, was introduced. With the assumption 
that the beam is inextensible, the strain enegy is expressed 


u - f </,.* ♦ (,) 

where * is the curvature of the beam. h and I 4 are geo 
metrical cross-sectional properties defined is follows. For 
edgewise deflection: 


h = 


12 

c 2l 

80 


and for flatwise deflection: 


( 2 ) 



where c and t are the height and thickness dimensions, re- 
spectively, of the cross section of the beam. Considering 
only the uncoupled, static deflections (load angles of 0 and 
90*), an equation for the deflection of the tip of the beam 
was derived as a function of the beam bending stiffness 
and other unknowns. Once the equation for the tip de- 
flection was derived, then a nonlinear least squares method 
was used to determine the best Eh and a that fit the ex- 
perimental data for both uncoupled flatwise and uncoupled 
edgewise deflections. The value obtained for a was ignored 
since GRASP does not consider material nonlinearity. The 
two values of E inferred from the bending stiffnesses and 
the cross-section geometry were averaged and multiplied by 
the cross-sectional area to obtain the axial stiffness. A value 
of Poisson’s ratio equal to 0.31 was assumed and the shear 
modulus, C, was inferred from E. The following stiffnesses 
resulted: 


axial stiffness = 6.2858 x 10* lb 
flatwise stiffness = 8.4487 x 10 3 Ib-in 3 
edgewise stiffness = 1.2689 x 10 4 lb-in 3 
torsional stiffness = 1.0538 x 10 3 lb-in 3 

The ratio of the edgewise stiffness to the flatwise stiff- 
ness should be 


77 = 15.97 (4) 

However, this ratio, based on the above stiffnesses reported 
above, turned out to be 


128.89 

8.4487 


= 15.02 


( 5 ) 


If on# act mu— that thia discrepancy la due to variations 
in width and thickness alone ths length of the beam, then 
the next question is bow much variation would it take to 
cause this discrepancy? Assuming an error t in each 
measurement, and substituting in the measured values for 
c and t, one obtains the ratio as 


(0.4090 + e) 2 
(0.1251 + 1) 7 


= 15.02 


(6) 


yielding t - 0.0052 in., a fairly small error. Thus, the 
inferred bending stiffnesses are not unreasonable. 

The moments of inertia of the various tip masses used 
were estimated with some gross assumptions since these 
values are relatively unimportant. The only properties of 
the tip masses stated in Ref. 9 and 10 were the masses. 
There was also a photograph depicting a tip mass with a 
hollow cylindrical shape. With this information in mind, 
several assumptions were made: the density of the tip mass 
was that of steel (0.284 lb/in. 3 ), the inner radius was 0.375 
in., and finally, the length of the tip mass was equal to its 
outer diameter. With these assumptions, the moments of 
inertia for tip masses of varying siie were calculated, (Table 
1 ). 


GRASP Model 


incident with the node TIP. Its definition specifies the mass 
and the mass momenta of inertias about all three principal 
axes. 


Correlation of GRASP Results With 
Experiment 

GRASP expresses static rotations in terms of Ro- 
drigues parameters , 12,3 so a min or amount of poetprocess- 
ing is needed to convert the GRASP output to the projected 
angle, (£), as measured in Ref. 0, ( p ^g. 3). Consider the or- 
thogonal triad at the root of the beam that remains aligned 
with the principal axes at the root. Introduce a dextral 
triad of unit vectors associated with those axes denoted by 
bf for i = l,2, and 3. Now consider a similar dextral triad 
at the tip of the beam, denoted by bf for « = 1,2, and 3, 
where the deflections and rotations were measured in the 
experiments. The relationship between the triads Is simply 

bf = C.ybf (7) 

where a repeated index implies summation. A line along the 
width of the croes-section is then aligned with bf = Ci,bf . 
Now consider the projection of bf in the plane determined 
by bf and bf denoted by p. The expression for p can be 
easily determined as 

P= b f-bf-b*b* 

=Cubf + Cub* 


The GRASP model for the Princeton experiment is 
depicted in Fig. 2 . Subsystem PRNCTN, the model-type 
subsystem generated internally by GRASP, represents the 
complete structure. The first explicitly defined subsystem 
is CANT BEAM. The frame of reference is defined to be co- 
incident with the model frame except for a rotation about 
the xj axis which is interpreted as the beam load angle. The 
subsystem cbntains two structural nodes named ROOT and 
TIP. ROOT ii coincident with the CANTBEAM subsys- 
tem frame of reference, and has all of its degrees of freedom 
prescribed to sero (cantilever beam boundary conditions). 
TIP io defined to be located 10.085 in. from the frame along 
the x 3 axis. 

The first child of CANTBEAM is an aeroelastic beam 
element named BEAM. An aeroelastie beam connectivity 
constraint associates the element 1 * root and tip nodes with 
the nodes ROOT and TIP in the subsystem CANTBEAM. 
The definition of the element includes specifying the or- 
ders of the polynomials used to represent the displacements. 
The typical approach in finite element programs would be 
to use several elements with the transverse displacements 
approximated by cubic polynomials, and the axial displace- 
ment and torsion approximated by linear polynomials. In- 
stead, for this analysis we use one element with eighth- 
order polynomials for bending and sixth-order polynomials 
for axial displacement and torsion. This yields a total of 
32 element degrees of freedom (8 of which are constrained 
out by the clamped-end condition). Essentially the same 
result! are obtained when the order of each polynomial is 
reduced by one. 

Subsystem WEIGHT, the second child of CANT- 
BEAM, is a rigid- body mas* element that is defined to be co- 


The angle measured in the experments is the angle between 
p and bf . From Fig. 3, it is clear that 


0 


.-•I 


Cx2 

Vl - Ci3 3 


In terms of Rodrigues parameters* 


( 9 ) 


C l2 




-*j + ^ 


( 10 ) 


where <£i, <£j, and <£3 are the Rodrigues parameters associ- 
ated with the rotation of the tip node. 

Also, all GRASP deflections have the deflections for 
no tip mass subtracted out before the results are plotted 
with the experimental data. All frequencies calculated by 
GRASP were converted from rad/sec to Hs. 

First, results are presented for a 1 -lb. tip mass. Figure 
4 shews the static deflections vs load angle. The GRASP 
correlation for flatwise and edgewise is excellent. Result* 
from Ref. 13 and 14 are shown here for comparison, with 
results from Ref. 14 presented only for the torsional de- 
flections. Transverse displacements from Ref. 14 were only 


* It should be noted that the matrix of direction cosines 
in thia work is the transpose of the one in Ref. 12 and that 
the Rodrigues parameters used in GRASP differ from those 
of Ref. 12 by a factor of 2 . 
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available for load angles of 30 and 40 # , and therefore, a 
complete load ai a sw e ep could not b« shown. It should 
be noted, however, that the transverse displacement results 
did agree well both with experiment and with GRASP. 

There were no dynamic results from Ref. 14. The cal- 
culations presented In Ref. 13 are based on the equations 
of Ref. 15. This analysis is restricted to moderate rotations 
caused by deformation, such that the squares of the rota- 
tional components are small compared to unity. Reference 
13 is not as accurate throughout the entire range of load 
angle as GRASP. For torsional deflection, the GRASP cal- 
culations cut right through the middle of the experimental 
scatter. The experimental scatter here is so large, however, 
that it is impossible to say which curve best fits the data. 
Figure 5 displays the flatwise and edgewise frequencies vs 
load angle. The GRASP results are only slightly offset from 
the experimental values, and follow the trend exactly. The 
average error is approximately 0.5%. Reference 13 does not 
pick up the trend for the flatwise frequency, however it does 
follow the trend for the edgewise frequency. Reference 14 
does not consider the dynamics. 

The 2-lb. tip-mass results are presented next. Here the 
torsional data have much lets scatter than in the 1-lb. case. 
Again, GRASP correlates excellently with the static deflec- 
tion as shown in Fig. 6. Also again, Ref. 13 is close but 
tends to deviate through certain portions of the load angle 
sweep. This deviation from the data is large for load an- 
gles above 40*. Reference 14, however, correlates quite well 
with the static data. Figure 7 shows the flatwise and edge- 
wise frequencies. The GRASP predictions are again slightly 
low for both of the frequencies, but follow the trends very 
nicely. Reference 13 while matching the data fairly well at 
0® load angle, strays from the data at higher load angles. 

Figure 8 presents the static results from the 3-Ib. tip- 
mass case. All three analyses appear to match the data 
well over the range shown. However, Ref. 13 results are 
available for only a load angle up to 15* , while Ref. 14 
calculated results up to 45*. The failure of Ref. 14 to 
converge past 15° is probably a result of the restriction to 
moderate rotation in that analysis. This is not a problem 
with GRASP, since GRASP has no such restriction. Ln 
Fig. 9 the flatwise and edgewise frequency sweeps are shown 
for the experiments! daU and the GRASP analysis. The 
correlation is excellent. 

For the 4-lb. tip-mass case there was not a very large 
experimental sweep available because of the large displace- 
ments which the beam underwent. GRASP did, however, 
correlate very well in the statics, Fig. 10, and the dynam- 
ics, Fig. 11. Reference 14 also correlates well with the static 
torsional deflection data. 

Figures 12 and 13 summarise flatwise and edgewise 
frequency vs tip mass. The lateral buckling load can be 
inferred from Fig. 12. As the tip mass increases, the flatwise 
frequency will tend toward xero. When lateral buckling is 
i mmin ent, the frequency will drop to xero very rapidly as 
indicated by both GRASP and Ref. 13. The edgewise 
frequency shown in Fig. 13 does not reach xero since there 
is no buckling in this mode. 


It should be noted that the analysis of Ref. 13 suf- 
fer* from being restricted to moderate rotations. Reference 
14 does much better than Ref. 13 in predicting the be- 
havior of this configuration because equations used therein 
are specialized for the type of structure used in the experi- 
ment. The equations of Ref. 14 are essentially identical to 
those of Ref. 13 except that certain terms of third degree 
in the unknowns are added to the analysis based on the 
observation that the coefficients of those terms are large. 
The size of those coefficients is a function of the ratio of 
the stiffnesses (thus depending on the cross-section geom- 
etry). These added terms would not be appropriate if the 
cross-section geometry were such that the stiffnesses were 
of the same order of magnitude. It is important to note 
that the equations in the GRASP analysis do not require 
that terms be added or removed in this msnner. This is 
an important consideration for general-purpose analyses, 
the equations for which should not have need of alteration 
merely because of changes in properties. 

Concluding Remarki 

GRASP is a general-purpose program with both the 
deUil and the generality to accurately model the end- 
loaded cantilever beam, as presented herein, ar well or bet- 
ter than the special-purpose analyses in Ref. 13 and 14. 
Although this experiment demonstrated significant nonlin- 
ear behavior both statically and dynamically, GRASP ac- 
curately predicts the results. The equations upon which 
GRASP is based are not restricted as far as the magnitudes 
of d£ j placement or rotation. Only the strains are required 
to be small compared to unity. The GRASP analysis is 
shown herein to be valid for these types of problems. As 
noted in the introduction, however, the present validation 
does not exercize many of the capabilities of the program. 
Also, as pointed out in Ref. 1, the analysis does need to be 
extended to treat Ixwms for which shear deformation would 
be important. 
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Table 1 Estimated Inertial Properties of Tip Mass 


Weight 

(»>•) 

Lateral Moments of Inertia 
(lb.-ln.-sec. 3 ) 

Axial Moment cf Inertia 
(Ib.-in.-sec. 3 ) 

1.0 

1.0822X10-* 

8.2356 x lO -4 

2.0 

3.3678x10-* 

2.6784 x 10~* 

3.0 

6.5673x10-* 

5.3169 x 10 _ * 

4.0 

1.0561x10-* 

8.6363 x 10“* 

5.0 

1.5276x10-* 

1.2573 x 10“ * 

6.0 

2.0658x10"* 

1.7083 x 10~* 

7.0 

2.6670x10“* 

2.2131 xl0 _ * 

8.0 

3.3278x10"* 

2.7691 xl0~* 

0.0 

4.0457x10“* 

3.3741 X10-* 

10.0 

4.8185x10-* 

4.0261 x 10“* 

10.42 

5.1589x10-’ 

4.3135 xlO"’ 

10.46 

5.1919x10-* 

4.3413 x 10“* 


i/i. t 
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Fig. 7 GRASP correlation with the Princeton experiment Fig. 8 GRASP correlation with the Princeton experiment 
(2-lb. tip mass): first flatwise and first edgewise frequencies (3-lb. tip mass): static deflections. 
y» load an*l«. .) fl.twU. v. load an 8 le. 

b) edgewise vs load angle. 

c) torsional vs load angle. 
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Fig. 0 GRASP condition with the Princeton experiment Fig. 10 GRASP correlation with the Princeton experiment 

(3-lb. tip mass): first flatwise and first edgewise frequencies (4-lb. tip mass): static deflections. 
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